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Although only a handful of organic materials have shown polariton condensation, their study is 
rapidly becoming more accessible. The spontaneous appearance of long-range spatial coherence is 
often recognized as a defining feature of such condensates. In this work, we study the emergence of 
spatial coherence in an organic microcavity and demonstrate a number of unique features stemming 
from the peculiarities of this material set. Despite its disordered nature, we find that correlations 
extend over the entire spot size and we measure g^\r,r') values of nearly unity at short distances 
and of 50% for points separated by nearly 10 /xm. We show that for large spots, strong shot to shot 
fluctuations emerge as varying phase gradients and defects, including the spontaneous formation 
of vortices. These are consistent with the presence of modulation instabilities. Furthermore, we 
find that measurements with flat-top spots are significantly influenced by disorder and can, in some 
cases, lead to the formation of mutually incoherent localized condensates. 

PACS numbers: 67.85.Hj, 81.05.Fb, 42.55.Sa, 71.36.+C 


In the last few years, a number of molecular sys¬ 
tems that show room-temperature polariton condensa¬ 
tion have emerged mu- Polaritons are hybrid exciton- 
photon quasiparticles that form in optical microcavities 
when dissipation is low as compared to the the light- 
matter interaction. Like their fundamental constituents, 
they obey Bose statistics at low densities. They inherit 
an effective mass as low as 10“ 10 times that of a Rb 87 
atom from their photonic component. Meanwhile, they 
collide with other polaritons and excitons due to an effec¬ 
tive interaction stemming from their matter component. 
In inorganic microcavities, these features have been ex¬ 
ploited to demonstrate rich phenomenology associated 
with Bose-Einstein condensation (BEC). The most well- 
known of these effects is the macroscopic occupation of 
the ground state beyond a critical density n c . A per¬ 
haps more important consequence of the phase transi¬ 
tion is the sudden appearance of off-diagonal long range 
order (ODLRO) in coordinate space [4j [5]. This mani¬ 
fests itself in the non-vanishing long-range behaviour of 
the first-order spatial coherence g^\ r, r’). At the po¬ 
lariton condensation threshold, a similar transition oc¬ 
curs, thus breaking U(l) symmetry, but also showing a 
number of distinct features resulting from the strongly 
non-equilibrium nature of the system [6]. 

Several materials have been used to demonstrate po¬ 
lariton condensation. The majority of these have used 
CdTe- and GaAs-based semiconductors whose operation 
is limited to low temperatures due to their small exciton 
binding energy. Recently, however, wide bandgap semi¬ 
conductors such as ZnO and GaN have emerged as vi¬ 
able materials for room-temperature applications [THE]. 
Organic semiconductors are especially attractive in this 
context due to their large exciton binding energy. Molec¬ 
ular excitons are commonly of the Frenkel type where 


both electron and hole are localized on a single molecule 
and exciton transport occurs incoherently between neigh¬ 
bouring molecules. In microcavities, long-range coupling 
between molecules is provided by the photon field and 
theoretical studies have shown that a phase transition 
can occur even in the absence of electronic coupling be¬ 
tween molecular dipoles m • Organics offer the advan¬ 
tage of a broad spectral range beyond that covered by 
GaN and ZnO and can be easily fabricated without the 
need for epitaxial growth. As-grown films, however, tend 
to be highly disordered and to overcome possible local¬ 
ization effects [T21I24] , the first demonstration of organic 
polariton condensation used a single-crystalline organic 
semiconductor as the active material pQ. Two recent 
demonstrations, however, have shown that the same phe¬ 
nomenology can be extended to amorphous systems: one 
consisting of a spin-coated polymer [31 and the other of 
a thermally evaporated oligomer [2]. By using interfero- 
grams, these reports have also evidenced the presence of 
long-range correlations in these systems. Here, we study 
in detail the onset of spatial and temporal coherence in 
an organic polariton condensate. We show that even in 
the presence of disorder, correlations can span the entire 
system size. Moreover, we highlight the important role of 
the pump and observe a number of unique features such 
as the spontaneous formation of vortices, which may sug¬ 
gest the presence of modulation instabilities. 

The microcavity is identical to that of Ref. [2]. It 
uses 9 dielectric mirror pairs on opposite sides of a 
layer of 2,7-bis[9,9-di(4-methylphenyl)-fluoren-2-yl]-9,9- 
di(4-methylphenyl)fluorene (TDAF) and was impulsively 
pumped high above the polariton energy (i.e. nonres- 
onantly). To highlight the presence of disorder, real- 
space photoluminescence (PL) using a Gaussian pump is 
shown in Fig. [l] In the linear regime, given in Fig. [lja), 
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FIG. 1. (color online), (a-d) Real-space images of the time- 
integrated PL below (a) and above threshold at different lo¬ 
cations (b-d) using an elliptical Gaussian pump. The linear 
image (a) remains identical regardless of the sample location, 
(e) Interferogram from location (d) averaged over 15 laser 
shots from, (f) Typical intensity versus phase delay. For ev¬ 
ery pixel, a cosine fit is used to extract both g^\— r,r) and 
the phase offset, (g) phase extracted from the fit. Scale bar 
is 5 /im. 


we find that the PL is independent of sample location. 
Above threshold, however, in-plane potential energy fluc¬ 
tuations cause strong variations as a function of sam¬ 
ple location, as shown in Fig. [ljb-d). Although elastic 
scattering due to disorder has been previously observed 
for organic polaritons in the linear regime [25] . the lo¬ 
cal intensity fluctuations seen here only appear beyond 
the condensation threshold. These bear striking similar¬ 
ities to those observed in disordered inorganic microcav¬ 
ities [IS El 123 GZ]. 

Although a number of other criteria exist, the spon¬ 
taneous onset of long-range spatial coherence is often 
recognized as a defining feature for condensation. This 
was first reported for a CdTe microcavity in the semi¬ 
nal work of Kasprzak et al using a Michelson interfer¬ 
ometer with one arm replaced by a retroreflector [27] . 
Soon after, the spatial coherence of a GaAs microcav¬ 
ity was measured using a Young’s double slit experiment 
and found to reach 30% at distances of 8 /im [28]. Theo¬ 
retical predictions regarding the asymptotic behaviour of 
g f 1 ) (r, r’) have shown that it can be used to shed light on 
the nature of the phase transition [29] and experiments 
with a low noise floor have observed a power-law decay 
and attributed this to a Berezinskii-Kosterlitz-Thouless 
(BKT) phase transition [50] . Recent theoretical work 
has, however, suggested that this decay may be the re¬ 
sult of a crossover phenomenon at intermediate length 
scales and that correlations should instead fall-off expo¬ 
nentially m- Interferometric techniques have also been 


used to identify excitations such as vortices and vortex- 
antivortex pairs and to probe quantum fluidic effects in 
polaritonic systems [321134] . 

In this work, we use the retroreflector configuration 
from Ref. a typical interferogram measured at 

the location of Fig. [ljd) is shown in Fig. [lje) . Despite 
the non-uniform intensity, parallel fringes indicate a flat 
phase over the entire condensate area. Similar results 
are observed for all sample locations. The visibility and 
phase offset are extracted by scanning a delay mirror 
(DM) over a phase of 67 r and fitting the resulting pat¬ 
tern at each pixel as shown in Fig. [ljf). The offset is 
shown in Fig. [ljg) and is due to the small angle differ¬ 
ence between the two interferometer arms. For a circu¬ 
lar spot, one expects ballistic propagation of polaritons 
away from the center accompanied by a radially increas¬ 
ing local wavevector for the condensate [55. An upper 
limit for this propagation length can be calculated us¬ 
ing l p < hk c /^ c m* = 2 /rm, where k c = 2.1 /1 m -1 is 
the largest condensate wavevector supported by the dis¬ 
persion, m* = 2.1 x 10 - 5 ra e is the effective mass and 
7 C = 6 ps -1 is the polariton decay rate. For small enough 
pump spots, this leads to an annular condensate density 
as shown in the Supplementary Material (SM). 

Figure [ 2 ] shows the spatial dependence of (— r, r) for 
3 values of pump fluence. The spatial coherence reaches 
peak values of g^\— r, r) > 0.8 for points separated by 
up to 5 /im. Vertical cuts, averaged over 8 pixels through 
the center of Fig. [2|a - c) are shown in Fig. §d - f). A 
vertical profile was chosen to avoid artefacts resulting 
from the small collapse along the horizontal direction. 
These occur due to the small difference in the arrival time 
of the 50° pump pulse along this axis. The noise remains 
significant near the condensate edge, which prevents a 
clear statement regarding the functional dependence of 
the asymptotic behaviour of g^\— r, r). In analogy with 
the thermal case, we simply show a Gaussian fit, which 
shows a reasonable agreement and allows us to extract 
a coherence length A c = 2\/2i rcr, where a is the stan¬ 
dard deviation of the Gaussian fit. We obtain a coher¬ 
ence length of up to A c = 38 /im, which is ultimately 
limited by the pump size. A Laical power dependence 
of g^(—y,y) is shown in Fig. |3[a). In all cases, an in¬ 
crease in 5 (1) ( —y,y) is observed, followed by a plateau 
near 1.5 P t h- 

As previously reported, the effective polariton- 
polariton and polariton-exciton nonlinearities in state- 
of-the-art organic microcavities are dominated by exci- 
ton saturation. The resulting change in refractive index 
causes anti-guiding in analogy with gain guided, index 
anti-guided conventional lasers. Moreover, the resulting 
nonlinearities are much weaker than in inorganics, which 
are dominated by the Coulomb exchange term between 
carriers [36]. From the blueshift in the region where the 
reservoir occupation is clamped (see SM), we can extract 
a polariton-polariton interaction coefficient g = 10 -6 
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FIG. 2. (color online), (a - c) Contour maps of g^\— r,r) for increasing pump fluence. The colorbar maximum corresponds 
to the peak value of g^\ r, —r) in the panel, (d - f) Vertical cuts through the (a - c) colour maps used to extract g^\—y, y). 
These are averaged over 8 horizontal pixels, with the standard deviation shown as error bars. The resulting coherence length 
is labelled A c and the fit residual is shown. Scale bars are 5 /xm. 


meV /im 2 . This allows us to estimate a typical healing 
length, calculated at 1.5 Pth, £ = ft/ y/m*g\^\ 2 = 2 /im, 
where |\k| 2 is the condensate density. This sets the scale 
for density fluctuations and is in reasonable agreement 
with our observations. Although the effect of disorder 
is apparent in the intensity patterns, the long coherence 
length and high degree of first-order coherence suggests 
the presence of a single condensate with a large conden¬ 
sate fraction [35] . 

For an equilibrium condensate, increasing density typ- 
a screening of the disorder potential. 
[22} have recently reported on the im¬ 
portant role of dissipation in non-equlibrium polariton 
condensates. They show that the competition between 
gain and dissipation can prevent screening of the disor¬ 
der expected when the polariton interaction dominates 
over the potential fluctuations. Similarly in this work, 
we do not observe a screening of potential fluctuations 
up to a blueshift of /x ^ 8 meV. 

The same interferometry setup was used to measure 
the first-order temporal coherence of the polariton con¬ 
densate. Here, the decrease in fringe contrast is mea¬ 
sured over long time delays for two points spatially sep¬ 
arated by 3 /xm. The temporal coherence, shown in 
Fig- §b), exhibits a Gaussian profile with a coherence 
time r c = 0.8 ps. Due to the impulsive pump, this 
value is not dominated by fluctuations, but instead agrees 
well with our calculation of the condensate survival time. 
From the emission linewidth, we find a condensate decay 
time of 1 ps, which is in good agreement with the mea¬ 
sured temporal coherence. 

One of the most startling features occurs for larger spot 
sizes. We find that for powers greater than ~ l.bPth 




FIG. 3. (a) The spatial coherence for increasing pump flu¬ 
ence. We find a sharp increase in —r), followed by a 

plateau near P = 1.5 Pth- Additional details are provided 
in the SM. (b) First-order temporal coherence, </d(r), for 
P = 1.7 Pth, measured for two points separated by 3 /xm. 
Values of g^\r) are shown as solid black circles and the er¬ 
ror bars indicate the standard deviation for each point. The 
temporal coherence of the condensate, r c , obtained from a 
Gaussian fit (red dashed line) is indicated as r c . 


phase gradients (density currents) and dislocations are 
commonly observed. These change from shot-to-shot 
and the overall interference pattern becomes washed out. 
This is shown in Fig. [4] where few-shot images were 
taken at the same location with small and large Gaussian 
pumps at l.SPth- Identical parallel fringes are observed 
for the small spot, but strong shot-to-shot variations oc¬ 
cur for the larger pump. For a small fraction of the laser 
shots, centro-symmetric fork dislocations due to singly 
charged vortices such as those shown in Fig. |4jc) can 
be observed. These form spontaneously at different lo¬ 
cations within the spot from shot to shot. Real-space 
profiles for a large pump also show density fluctuations, 


ically leads to 
Thunert et al. 
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FIG. 4. (color online), (a - b) Interferograms using a small 
(27 x 25 jttm) Gaussian pump taken at the same sample loca¬ 
tion and power integrated over 4 laser shots, (c-d) Interfero¬ 
grams using a large (97 x 60 /xm) Gaussian pump integrated 
over 2 laser shots. Yellow circles indicate a fork dislocation. 
Scale bar is 5 /xm. 



FIG. 5. (color online). Real-space emission for a flat-top 
pump spot below (a) and above threshold (b - d). (c) and 
(d) are intensity patterns that were recorded from different 
sample locations, (e - f) interference patterns recorder at the 
same pump power (1.5Pt^) and represent the sample locations 
(b) and (c), respectively. Note the change in fringe orienta¬ 
tion, cross-fringes and poor contrast. Scale bar is 15 /xm. 


while occasional regions of vanishing density can be at¬ 
tributed to vortex cores (see SM). 

To emphasize the crucial role of the pump shape, we 
have performed similar experiments using a refractive 
beam shaper to control the pump beam profile. In this 
case, the excitation spot was made significantly larger 
and with a flat-top intensity. Figure [5] shows the below 
threshold (a) and above threshold (b - d) real-space im¬ 
ages from 3 different sample locations. Below threshold, 
the intensity is nearly flat over the entire pumped region, 
albeit with some small fluctuations. Above threshold, 
the appearance of multiple bright regions is observed. In 
some cases, these regions are separated by as much a ten 


microns. The resulting interferograms tend to be highly 
distorted often showing a mixture of dislocations, phase 
gradients and incoherent regions with strong shot-to-shot 
fluctuations. We have, in addition, verified that in cases 
with a significant separation between condensates, mu¬ 
tual coherence is lost (see SM). 

Our experiments with various pump sizes and shapes 
point to two distinct causes for these distorted patterns. 
First, a greater sampling of the disorder for large spots 
can lead to condensation in regions further away from 
the center. In this case, our interferometric geometry 
which inverts about the origin will sample regions that 
can be uncorrelated or possess additional phase shifts 
due to their different positions. The latter is seen in 
Fig. [5je-f) where the fringe orientation changes despite 
the identical geometry. Note also the presence of smaller 
cross-fringes in both cases. These features are never seen 
for Gaussian pumping. Second, small Gaussian pump¬ 
ing has been theoretically predicted to be beneficial in 
reducing the region in reciprocal space where modula¬ 
tion instabilities can occur m [37H39]. The criteria for 
such instabilities to occur is nearly always satisfied in 
organic cavities, which possess typical cavity decay rate 
> 10 2 times those of reservoir excitons. In our case, for a 
homogeneous spot, these instabilities exist over a broad 
range of wavevectors for pump powers up to 10 4 Pth- The 
short survival time of our condensate helps prevent their 
significant growth. 

In conclusion, we find that organic polariton conden¬ 
sates bear similarities to disordered inorganic microcav¬ 
ities, but with some distinct features. For smaller spots, 
we always observe a single condensate showing good 
phase correlation and a coherence length ultimately lim¬ 
ited by the pump size. In contrast, for large Gaussian 
pumps, phase gradients and defects emerge, with strong 
shot-to-shot fluctuations. Although masked in the linear 
regime, the role of disorder becomes evident beyond the 
condensation threshold. For homogeneous spots, the in¬ 
terferograms become highly distorted as a result of both 
disorder and fluctuations. The size and profile depen¬ 
dence are consistent with current predictions and suggest 
that the fluctuations may be due to modulation instabil¬ 
ities. Another interesting possibility, however, is that 
these may be remnants from a Kibble-Zurek-like mech¬ 
anism [39]. The ease of fabrication of TDAF — it can be 
spin-coated or thermally evaporated — make it a useful 
platform for the further investigations of these effects. 
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SUPPLEMENTARY MATERIAL 
1. Sample and Experimental configuration 

As in our previous work, the microcavity used in this Letter consists of a 120 nm thick TDAF layer 
sandwiched between two 9-pair distributed Bragg reflectors (DBRs). These consist of alternating 
SiC >2 and Ta 2 Os layers (Fig. Sir). The sample is excited non-resonantly using 250 fs pulses from an 
optical parametric amplifier tuned to 3.22 eV. The pulses are TM-polarized and incident at 50° with 
a repetition rate of 1 kHz for PL and time-integrated interferometry and 100 Hz for the single-shot 
interferograms. Figure [Slj) shows the TM-polarized PL measured on a location with a detuning A 
= -400 meV. The dispersion parameters are indicated in the figure. 


(a) 

9x Si02/ 

Ta 2 0 5 

TDAF 

9 ** 020 - 

Ta 2 0 5 



0 10 20 30 40 

Angle (deg) 


FIG. SI. (color online), (a) Schematic of our microcavity, (b) TM-polarized PL of a negative detuned cavity (A = -400 meV). 
The detuning was obtained using the coupled harmonic oscillator model described in Ref. [2J. Lower polariton (LP) and cavity 
mode (C) dispersions are shown with dashed lines and the fit parameters are indicated 


Our experimental configuration is shown schematically in Fig. |S2|(a). In the retroreflector Michel- 
son interferometer configuration, when the image centers are aligned, the first-order spatial coherence 
between centrosymmetric points g^ l \ r, r’) = g^\ r, — r) can be obtained from the visibility, V(r), of 
the interferogram fringes: 


i»= 


1-max I min 
1-max I min 


2y/(r)/(—r) (1) 


/(r) + /(-r) 


9 ( TU-r) 


Where /(r) is the light intensity and I ma x and I mm are respectively the maximum and minimum 
intensities measured, when the phase delay is varied. 


The power dependence and blue shift measured for our microcavity are given in Fig. S2(b). The 
threshold energy is P t h = 600 ±100 pJ, which corresponds to an absorbed pump fluence of 27 ± 
5 /.(J/cm 2 . The errors correspond to the standard deviation over various sample locations. 


2. Large-scale collapse 

Beyond threshold, we observe a collapse of the emission region towards the pump center (~ 70% 
reduction in full width at half maximum). A typical vertical intensity profile is shown in Fig.|S2[c) as 
a function of increasing pump power. We find that this large-scale collapse can be readily explained 
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FIG. S2. (color online), (a) Schematic of our experiment where one arm of a Michelson interferometer is replaced by a 
retroreflector (RM). A piezoelectric translation stage is to vary the phase introduced by the delay mirror (DM). The reflected 
images from RM and DM are then overlapped on a CCD camera (Thorlabs BC106N-VIS). (b) Power dependence of the PL 
and blueshift (inset) observed at normal incidence, (c) A typical vertical profile of the intensity for a small Gaussian pump 
as a function of increasing pump fluence. (d) Profile calculated by scaling the linear regime profile by the peak PL power 
dependence. The inset shows the profiles on a linear scale. 


by scaling the intensity profile in the linear regime, which follows the pump profile, by the normal- 
incidence power dependence measured using a fiber-coupled CCD spectrometer shown in Fig. S2(b). 
If n(r) is the normalized pump profile in the linear regime, f(p ) the power dependence at power 
p = P/Pth, the intensity is calculated such that /(r) = f(n(r)p). The agreement obtained is 
illustrated in Figure |S2fd). Note, however, that spatial disorder is not captured by this model. The 
non-uniform density peaks within the condensate apparent in Fig. |S2fc) are absent in the linear 
profile and consequently in Fig. S2(d). 


3. Outwards propagation of polaritons 

For small enough pump spots, the above threshold emission becomes distorted into a ring-like 
pattern in contrast to the homogeneous large-scale collapse. We attribute this to the outwards flow 
of polaritons due to the repulsive potential induced by the pump. This behaviour is shown in Fig. [S3] 
and qualitatively observed everywhere on the sample. Disorder typically causes some deviation from 
a perfect ring-shape and we have carefully verified that this behaviour is reversible and not a result 
of damage within the central region. 


4. Power-dependence of spatial coherence 


We have measured the power-dependence of gP(—y, y ) on three different sample locations using 
a small Gaussian pump. These are shown in Fig. S4 around gP( 0,0) and at distances \2y\ = 4 //in 
and \2y\ = 8 /on. In all cases, values and standard deviations were obtained by averaging over 25 
pixels on both sides of the origin. 









































FIG. S3, (color online). Typical real-space condensate emission measured as a function of increasing pump power for a small 
Gaussian pump 



FIG. S4. (color online), (a - c) The spatial coherence at various locations away from the condensate center r = 0 fim. shown 
for increasing pump fluence. Note the sharp increase in g^\—y, y), followed by a plateau near P = 1 .5Pth, The value and error 
bars correspond to the mean and standard deviation obtained from averaging over 25 pixels. 


5. Additional real-space and single-shot images 

Figure [S5] shows supplementary few-shot interferograms images measured for the three different 
spot types. Averaged real-space images at the same location are included. Note the inhomogeneous 
character of the large Gaussian real-space image. We can also see that the flat-top pump results in a 
significantly larger number of phase defects as compared to a Gaussian pump. Single-shot real-space 
images comparing the behaviour for small and large pumps are shown in Figure [S6| Although we 
cannot simultaneously acquire a single-shot real-space image and interferogram in our setup, the 
occasional appearance of zero intensity regions in the real-space image with dimensions on the order 
of the healing length is consistent with the vanishing density of a vortex core. 


6. Mutual coherence between independent condensates 

In some cases, flat-top pumping leads to regions of high density separated by as much a ten microns. 
To probe whether these various regions are mutually coherent, we investigated a sample location, 
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FIG. S5. (color online), (a-c) Real-space (a) and 4-shot interferograms (b-c) for a small Gaussian pump, (d-f) Real-space 
(d) and 2-shot interferograms (e-f) for a large Gaussian pump, (g—i) Real-space (g) and 2-shot interferograms for a flat-top 
pump (h-i). Real space images were integrated over a) 100, d) 30 and g) 80 laser shots. Scale bar is 5 /xm. 



FIG. S6. (color online). Two consecutive single-shot images taken above threshold for (a-b) a small Gaussian pump and (c-d) 
a large Gaussian pump. Scale bar is 5 /xm. 


given in Fig. |S7ja), showing 3 seemingly isolated condensates denoted as 1, 2 and 3. By translating 
the mirror and retro-reflector images, we have performed interferometric experiments overlapping 
the central spot 2 with itself (b), with spot 1 (c) and spot 3 (d). Clear parallel fringes are observed 
in Fig. |S7fb) indicating that this condensate still shows strong phase correlations. However, no 
fringes can be resolved in Fig. |S7fb) and (c), which correspond to the interference between different 
condensates. This evidences the presence of three distinct, but mutually incoherent condensates. 
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FIG. S7. (color online), (a) Real-space image for a large pump showing 3 localized condensates denoted as 1, 2, 3. (b - d) 
Interference pattern obtained by manually overlapping condensate 2 with itself (b), with 1 (c) and 3 (d). Scale bar is 15 /im. 






